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ABSTRACT 


This study examines the integration of ocean glider Unmanned Surface Vehicles 
(USVs) and Unmanned Undersea Vehicles (UUVs) in support of wide-area 
oceanographic and acoustic sampling. These collaborative systems could enable the U.S. 
Navy to conduct multi-month unmanned maritime surveillance. Optimal sensor position 
in the water column and persistence are critical requirements to reduce surface 
expressions of such a network. An experiment was conducted in Monterey Bay to 
evaluate underwater gliders as mobile passive acoustic sensing platforms. Acoustic 
propagation modeling was used to plan experiment geometry, predict transmission loss 
(TL), and estimate acoustic communications performance with a USV. A medium 
frequency acoustic source was deployed at a range of 5.5 km from a receiver on board a 
glider conducting a 1000 m dive to demonstrate that a glider can adapt to the local 
environment to exploit more reliable propagation paths. Results demonstrate that gliders 
are effective mobile platforms to support persistent acoustic sensing. The glider received 
transmitted signals at levels in close agreement with TL predictions. Signals were 
received while the glider was in motion, and reception improved during a quiet deep 
loiter. Given the depths, ranges, and environmental conditions studied, research and 
modeling suggest sufficient acoustic communication performance to promote 
connectivity of the proposed network. 
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I. 


INTRODUCTION 


The long-endurance of the autonomous underwater glider platform makes it an 
excellent candidate to support a mobile undersea acoustic cueing system. Such systems 
could enable the U.S. Navy to conduct persistent autonomous battlespace surveillance 
and reconnaissance. While this concept might seem far-fetched, it is supported by 
numerous technologies that are currently in use with underwater glider research. Some of 
these advancements include compact acoustic towed and tetrahedral arrays (Hughes et al. 
2009, Maguer et al. 2013), onboard signal processing (Dassatti et al. 2011), and 
automated acoustic signal classification (Baumgartner et al. 2014). 

This thesis evaluates underwater gliders as mobile, autonomous, passive acoustic 
sensing assets. Of specific interest are the ability of these platforms to loiter at depth and 
to exploit the reliable acoustic propagation paths present in a deep ocean environment. In 
this study an autonomous underwater glider functions as a node in a theoretical 
collaborative network with an unmanned surface vehicle (USV) that operates as a 
wireless acoustic communications relay, or gateway, node. Research demonstrates this is 
a functional network model with potential for future use (Bingham et al. 2012). 

A thorough analysis of how to design and operate such a system is warranted. 
This study, however, considers two critical aspects of the inherent challenges: where to 
position the sensing platform in the water column and how to keep it there for an 
extended loiter. A field experiment and acoustic propagation modeling were performed to 
investigate these operational factors. Data analysis and results are provided to evaluate 
acoustic sensing performance during phases of a glider dive cycle and a deep loiter. 
Modeling of acoustic channel impulse response is included to support a theoretical 
analysis of acoustic communication performance for the proposed system. 
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II. BACKGROUND 


The U.S. Navy’s Information Dominance Roadmap for 2013-2028 predicts 
innovative and covert solutions are needed to effectively counter modem adversaries that 
employ anti-access/area-denial (A2/AD) strategies (U.S. Navy 2013). One suggested 
solution is to develop autonomous vehicles that can gain access to denied environments 
and conduct environmental sensing. The report states that success in the “information¬ 
intensive combat environments of the future” will in part depend on the capacity of these 
vehicles to enhance battlespace awareness. The roadmap defines Battlespace Awareness 
(BA) as, “the ability to understand the disposition and intentions of potential adversaries 
as well as the characteristics and conditions of the operational environment.” The report 
also identifies major elements of BA that the Navy needs to optimize in order to 
effectively meet future threats. Requirements for improved data collection and processing 
stand out from that list as elements of BA that will benefit from the contributions of 
unmanned systems. Autonomous vehicles that operate in the undersea domain are 
expected to play a vital role in this collaborative effort. 

A. MOTIVATION FOR GLIDER USE 

It is essential to select the ideal vehicle design to meet specific undersea mission 
requirements. The 2004 edition of the U.S. Navy’s Unmanned Underwater Vehicle 
(UUV) Master Plan identifies vehicle classes that are deemed appropriate for a host of 
proposed missions. It recommends underwater gliders for the collection of environmental 
data in regions, “where battlespace dominance has not been achieved” (U.S. Navy 2004). 
Underwater gliders are a type of unmanned undersea vehicle (UUV) that use a buoyancy 
engine to move vertically and wings to move horizontally throughout the water column 
(Rudnick et al. 2004, Sherman et al. 2001). Underwater gliders use their sensors to 
characterize the oceanographic environment and, due to their unique method of 
operation, have far greater range and endurance than other unmanned undersea vehicles 
that utilize more conventional propulsion (Rogers et al. 2004). Proven designs, such as 
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the Slocum glider (see Figure 1) and the Spray glider (see Figure 2), are two-man 
portable as well as simple to deploy and recover. 

Figure 1. The Slocum glider developed by Teledyne Webb Research as it is 

launched from a small boat. 



Slocum gliders have a modular payload design (Teledyne Webb Research 2015). Various 
models can operate at depths from 4-1,000 meters, for up to 365 days, with a max range 
of 6,000 km (from Naval Drones 2015, http://www.navaldrones.com/images/glider.jpg). 


The UUV Master Plan also suggests that pre-positioned autonomous cueing assets 
will support future Anti-Submarine Warfare (ASW) and Intelligence, Surveillance, and 
Reconnaissance (ISR) missions (U.S. Navy 2004). The use of gliders for acoustic signal 
detection and classification is not specifically mentioned. However, in the decade since 
the publication of the 2004 UUV Master Plan these capabilities have been consistently 
demonstrated with various glider platforms. When outfitted with underwater acoustic 
receivers, known as hydrophones, and onboard signal processing capabilities gliders 
become quite effective mobile passive acoustic sensing platforms (Baumgartner et al. 
2014). Regardless of the specific type of sensors, whether acoustic or non-acoustic, it is 
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quite easy to imagine how networks of such vehicles will help facilitate the transition 
from legacy methods of undersea environmental assessment. 


Figure 2. The Spray glider. 



Cutaway diagram of the Spray glider developed by the Scripps Institution of 
Oceanography (SIO) and the Woods Hole Oceanographic Institution (WHOI). The 
battery packs, inside a sealed pressure hull, are mechanically shifted to control vehicle 
pitch and roll attitude. The external bladder, within the flooded science bay, fills 
with oil to change vehicle buoyancy. Antennas for satellite communications and GPS 
navigation, not depicted, are built into the glider’s wings. Mission endurance can last six 
months and cover 3,600 km (from Scripps Institution of Oceanography 2015, 
http://spray.ucsd.edu/pub/rel/info/spray_description.php). 


B. N2/N6 PROPOSED STUDY 

Given the capabilities of currently available platforms, and the potential to 

execute missions with military relevance, a study of the integration of US Vs and UUVs 

in support of wide-area environmental sensing and characterization has been proposed by 

the Deputy Chief of Naval Operations for Information Dominance (DCNO N2/N6). In 

support of this study a Naval Postgraduate School research proposal was submitted that 

outlines the concept of joining USV and UUV assets into a mobile, collaborative, and 

long-endurance network for both oceanographic and acoustic environmental 

characterization (Joseph and Horner 2014). The authors specifically hypothesize that, 

“the most critical component for effective UUV/USV integration is reliable cross-domain 

communications to effectively exchange navigation, command and control (C2) and 
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environmental information.” This thesis expands on research questions presented in the 
proposal that encourage an evaluation of limitations that could impact the effective 
performance of such a system. The theoretical network framework described in the 
proposal is the foundation for the discussions presented in this study. 

C. CONCEPT OF OPERATIONS 

Conventional glider operations require routine trips to the surface to transmit 
recorded sensor data and vehicle status reports. Once on the surface, gliders establish a 
satellite communication link and report the requested data to distant controlling centers. 
Routine trips to the surface are also required for GPS signal reception in order to receive 
position updates and maintain accurate navigation. If a network of underwater gliders 
with passive acoustic sensing capabilities were to remain submerged, and thereby reduce 
the surface expressions of such a system, the vehicles would need a reliable and efficient 
method to transmit crucial data and receive position updates. This problem could 
potentially be solved through collaboration with a GPS-enabled USV and acoustic 
communications. This would allow an autonomous asset to persist at an ideal depth for its 
select suite of onboard sensors. 

1. OPTIMAL SENSOR PLACEMENT 

As acoustic energy propagates away from an underwater sound source it follows 
many paths that are individually influenced by physical parameters that change with 
range and depth. Variations in parameters such as temperature, salinity, and pressure 
affect how fast sound travels through seawater. As propagating sound waves encounter 
regions of the ocean with different sound speeds the waves respond by refracting. In the 
mid-latitudes and tropics the strength of this refraction is driven by the relatively rapid 
decrease in water temperature with depth. This occurs in a layer of the ocean known as 
the main thermocline that extends below the warmer surface layer. As ocean depth 
increases in this layer the rate of change for sound speed, known as the sound speed 
gradient, is negative and sound is refracted downward as it propagates away from a 
source. When a portion of the total radiated acoustic energy from a sound source follows 
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this downward refracting path to a receiver, without first interacting with the ocean 
surface or bottom, it is propagating along a Reliable Acoustic Path (RAP) (see Figure 3). 

Figure 3. Reliable Acoustic Path. 



Depictions of various reliable acoustic paths that can exist between a deep acoustic 
source and a near-surface receiver (from Urick, R. J., 1983: Principles of Underwater 
Sound: Third Edition. Peninsula Publishing, 423 pp.) 

The path is referred to as reliable because it is a direct path, from source to 
receiver, that avoids the boundary interactions that redirect and degrade acoustic signal 
intensity (Urick 1983). If a homogenous ocean is assumed then acoustic reciprocity 
implies this path can be exploited in reverse by simply alternating the locations of 
the source and receiver (Rayleigh 1894). This makes the RAP a valuable physical 
phenomenon of the deep ocean for both passive acoustic detection and acoustic 
communications. Given the low ambient noise levels at greater depths this enables a deep 
acoustic sensor to monitor a much larger area at a high signal to noise ratio (SNR) 
(Thompson 2009). 

Because of these factors, there is great utility in having a sensing platform that 
can position itself at an ideal location in the water column, engage in a multi-month 
loiter, and sample the environment while it gets there. Optimizing sensor placement in 
this way increases acoustic detection ranges and time in contact. The field experiment 
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conducted for this study seeks to demonstrate these benefits by conducting environmental 
sampling in transit and during quiet loitering. 

2. LONG-ENDURANCE LOITERING 

Gliders and other undersea vehicles modified with underwater acoustic telemetry 
modems have demonstrated that trips to the surface to transmit data can at least be 
minimized when a USV (see Figure 4) acts as a network surface gateway platform 
(Maguer et al. 2013, McGillivary et al. 2012). Similar research has also shown that 
acoustic communications make it possible for a USV to send commands (Howe et al. 
2010), as well as estimate the range and bearing to underwater gliders and thus serve as a 
navigation reference gateway (Bingham et al. 2012, Stanway et al. 2014). 

Therefore, underwater acoustic communications provides one method to connect 
a distributed network of autonomous undersea and surface vehicles and enable long- 
endurance loitering. According to this system model, the portion of a glider’s battery 
power normally allocated to operate a buoyancy engine would be available for tasks such 
as onboard signal processing and acoustic communications. This would be significant 
considering that it is not uncommon for conventional glider designs to allocate more 
than eighty percent of total battery power to buoyancy engine operation (Jenkins et al. 
2003). A glider with passive acoustic sensing capabilities could then loiter at an optimal 
depth determined by its own environmental sensing, process data onboard, and quickly 
communicate critical information off-board without the need for trips to the surface. A 
distributed network of these platforms would then have a limited surface expression; one 
USV could coordinate with many gliders. Such a network has great potential to become a 
valuable asset for enhancing battlespace awareness. 

While this method of collaboration can potentially reduce the long 
communication time-delays imposed by glider dive cycles, underwater acoustic 
communications is not a perfect solution. The underwater environment between any 
acoustic source and receiver, referred to as the acoustic channel, places limits on 
available transmission bandwidth. Transmitted signals must also contend with various 
physical challenges of the time varying undersea environment such as surface conditions, 
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bathymetry, and ambient noise levels (Rogers et al. 2004, Howe et al. 2010). As a result 
the rate that data successfully arrives at a receiver modem, known as network throughput, 
is much lower than more traditional methods for connecting communications networks. 
Therefore, acoustic communications must be utilized in an efficient manner in order to be 
a viable method of connecting such a network. 


Figure 4. Autonomous Surface Vehicles (USVs). 



(a) (b) 


Examples of autonomous unmanned surface vehicles (USVs) built by Liquid Robotics: 
(a) exploded view of the hybrid Wave Glider SV3 (from Liquid Robotics 2015, 
http://liquidr.com/technology/waveglider/sv3.html); (b) silhouette that depicts a Wave 
Glider SV2 modified to operate an acoustic modem housed in a towed body (from Liquid 
Robotics 2015, http://hquidr.com/resdown/resources/case-studies/cornell.html). The 
vehicles are solar-powered and utilize a customizable surface float to carry an assortment 
of computer payloads for environmental sensing, communication, and navigation. A 
tethered sub-surface platform with wings harnesses wave energy to propel the glider to 
waypoints. 
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III. OCEAN ACOUSTICS 


In order to understand why autonomous vehicles can use acoustics to passively 
sense the underwater environment and communicate with acoustic modems it is 
important to first understand how sound moves through water. Sound is a mechanical 
longitudinal wave of energy. A sound wave is set in motion in an elastic medium, such as 
seawater, by a vibrating source that creates pressure disturbances. In other words, sound 
travels through the ocean as pressure waves. Sound propagation in seawater can be 
described mathematically through solutions to the acoustic wave Equation (1) using 
appropriate descriptions of the boundaries and medium (Urick 1983). 


c 2 dt 2 


( 1 ) 


The acoustic wave equation describes a three-dimensional change in acoustic 
pressure p with respect to position, time /, and sound speed c . Solutions to the equation 
explain how much of the original sound pressure remains after a signal travels a certain 
distance from a source. This technique enables computer models to numerically predict 
the motion and decay of acoustic wave energy in the ocean. 

A common theoretical approach used to solve the acoustic wave equation in these 
models is known as ray theory. This geometric approach uses rays to trace the spatial 
path, or trajectory, of radiated acoustic pressure. The amplitude of acoustic pressure 
along the trajectory of a particular ray is then represented by the cross-section of a ray 
tube formed by adjacent rays. A detailed derivation of the ray theory approach and its 
application in computer-based models can be found in Jensen et al. (1994). A general 
overview of the theory and how it models sound propagation is provided below. 


A. RAY THEORY 

As previously described, the speed of sound in seawater plays a central role in 
determining how acoustic energy propagates. Sound speed in the ocean is a function of 
temperature, salinity, and pressure. Representations of its variation with respect to depth 
are called sound speed profiles. These variations in sound speed cause rays of acoustic 
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pressure energy to follow curved paths from source to receiver. This behavior is known 
as refraction and is described by an acoustic form of Snell’s Law (Hovem 2013) 

^ = cos6>(z) = cosfl 0 (2) 

c(z) c 0 


While sound speed depends on depth c(z), it has range-dependence as well. 
However, it is more straightforward to illustrate Snell’s Law from a range independent 
perspective. If a medium is divided into horizontal layers with constant velocity the law 
explains that acoustic ray grazing angles at layer boundaries are directly influenced by 
the speed of sound in those layers (Urick 1983). As a result, incident rays change 
direction upon crossing a layer boundary. As these direction changes accumulate over 
distance and time, the ray path takes on the shape of an arc. This can be visualized by 
considering a particular segment of an arcing acoustic ray path (see Figure 5). The shape 
of the ray segment is defined by a change in range and depth, dr and dz , and its angle 0 
relative to the horizontal. When sound speed varies with depth, known as the sound speed 
gradient g(z), 




dc{z) 

dz 


(3) 


the ray angle 6 will change according to Snell’s Law. Ultimately, this means that the 
radius of curvature R(z) for this ray path. 


*(z) 


c ( z ) 1 

COS e(z) g(z) 


1 

Zs(z)’ 


(4) 


is also a function of depth (Hovem 2013). The ray parameter E, defined in Equation (2) is 
a constant ratio. It is a unique value for a particular ray. The curvature of an acoustic ray 
is simply a function of this constant value and the sound speed gradient at a specific 
depth. The sign of the sound speed gradient indicates if the ray will curve upward or 
downward. This means that if sound speeds are trending slower with increasing depth, a 
negative gradient, acoustic rays will refract downward (see Figure 6). Ray tracing 
computer programs use these techniques to incrementally map out individual ray paths 
and calculate acoustic travel times to receiver locations. 
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Figure 5. Arcing acoustic ray path 


r 



A ray path segment with an arc length ds in an isotropic medium (from Hovem, J. M., 
2013: Ray trace modeling of underwater sound propagation. Modeling and Measurement 
Methods for Acoustic Waves and for Acoustic Microdevices, M. G. Beghi, Ed., InTech, 
DOI: 10.5772/55935). 


Figure 6. Diagram of an acoustic ray trace 




(a) (b) 

(a) a negative gradient sound speed profile; (b) representations of how acoustic rays with 
various initial angles will propagate from a source depth Z t given the profile, as well the 

depiction of a “shadow zone” that forms at the limit of the tangential ray path (from 
Brekhovskikh, L. M., and Yu. P. Lysanov, 2003: Fundamentals of Ocean Acoustics. 
Springer, 280 pp.). 
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Once an ensemble of acoustic ray paths are calculated, ray tracing models must 
also calculate how acoustic pressures will change along the way. This is modeled by the 
varying cross section of a tube formed by adjacent acoustic rays (see Figure 7). This ray 
tube is considered to have constant acoustic power and its width dL is then proportional 
to the amplitude of acoustic energy at some distance from a sound source. The average 
acoustic pressure, or root mean square pressure p rms , that moves through a unit cross- 

section of the tube for a unit time represents acoustic power flux and is called acoustic 
intensity I 

2 

I _Prms_ , (5) 

pc 

where p and c are seawater density and sound speed, respectively. 


The “loudness” of sound at a location is simply a ratio of two intensities or 
pressures. When the ratio utilizes a standard reference intensity or reference pressure the 
result is a parameter known as sound pressure level SPL 


SPL = 10 log 
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where I ref =6.7x10 — j, and p ref =\pPa . This “loudness” is represented by the 


m 


dimensionless logarithmic unit decibel ( dB ) and the associated reference parameter. 


When sound pressure levels are calculated in the ray tube at one meter from a 
sound source, as well as at some distance r from the source, the comparison yields the 
estimated reduction in sound level due to propagation, or transmission loss TL 

(7) 
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TL = 20 log 
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However, this is only a geometric representation of how sound pressure changes 
as it propagates from a source. There are several other important factors to consider that 
influence transmission loss and they are commonly implemented in ray theory models. 
These will be discussed in the section that addresses the passive sonar equation. 
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Figure 7. Ray tube model 



Acoustic intensity modeled with ray tubes. Range and depth increase along the horizontal 
axis and vertical axis, respectively. The gray circle represents a sound source with initial 
acoustic ray angle 6 {) ; cl6 0 is the initial separation between two rays; dr is an 

increment in range; dz is the depth difference between the two rays; and dL is the ray 
tube width (from Hovem, J. M., 2013: Ray trace modeling of underwater sound 
propagation. Modeling and Measurement Methods for Acoustic Waves and for Acoustic 
Microdevices, M. G. Beghi, Ed., InTech, DOI: 10.5772/55935). 


B. MODIFIED RAY TRACING 

While the traditional ray tracing method is computationally fast, it has several 
important limitations. It produces excellent results for higher acoustic frequencies 
(e.g., > 1 kHz), but fails at lower frequencies. This is due to the fact that ray tracing 
produces approximate solutions to the wave equation. It does this through an assumption 
that the scale of the sound speed gradient is much larger than the acoustic wavelengths. 
This is fine for higher frequencies where wavelengths are short. Lower frequencies, 
however, with longer wavelengths make this critical assumption inappropriate (Urick 
1983). Other models based on normal mode theory and the parabolic equation method 
produce exact solutions for a wider range of frequencies but can be computationally 
expensive. Approximations produced by ray tracing also incorrectly predict infinite 
acoustic intensity or no acoustic intensity at regions known as caustics and shadow zones, 
respectively. However, a modified ray tracing model known as Bellhop overcomes these 
limitations and produces predictions in excellent agreement with exact solution models 
(Porter 2011). This model was chosen for use in this study, but it is important to note that 
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it is not a complete replacement for exact solution methods. Its strengths and limitations 
are explained in the next chapter. 

C. THE PASSIVE SONAR EQUATION 

Acoustic sensing and acoustic communications are fundamentally similar in 
that both rely on the successful detection of sound at some distance from a source. Ray 
theory shows that acoustic signal energy will dissipate as it propagates along various 
underwater paths. Acoustic receivers must then detect the presence of this signal amid the 
background noise. Detection is therefore an analysis of the input signal-to-noise ratio to 
uncover signals that sufficiently exceed the noise level (Urick 1983). This process 
requires an application of a threshold that once exceeded results in a conclusion that the 
signal of interest is present in the noise. The threshold must strike an important balance 
between two competing factors: probability of detection and probability of false alarm. 
This is a challenge because when the threshold is too high both the probability of 
detection and the probability of false alarm are low. When set too low, both probabilities 
are high (Urick 1983). 

The proper balance depends on the environment, sonar system design, and signals 
of interest. Signals and noise, however, are affected by the environmental variability in 
the ocean that induces acoustic uncertainty. As a result, these vital parameters constantly 
fluctuate from one instant to the next. Each parameter is a stochastic, or random, process 
that can only be estimated. 

In order to predict the performance of an acoustic detection system it is necessary 
to establish a relationship between the passive sonar parameters: source level ( SL ), 
transmission loss (TL ), noise level ( NL ), directivity index ( DI ), and detection threshold 
(DT ). This working relationship is known as the passive sonar equation (Urick 1983). 

SL-TL- NL-DI + DT (8) 

It is important to keep in mind that this simple equation represents a relationship 
among parameters with inherent uncertainty and does not reduce underwater acoustics to 
simple deterministic solutions. The directivity index (DI) describes the ability of the 
receiver to resolve acoustic detections from a particular bearing. In the case of this study, 
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an omnidirectional receiver is used so DI can be set to zero. The detection threshold 
(DT ) represents the minimum remaining signal at a receiver, in excess of the noise level 
( NL ), that is required for a detection. This implies that the passive sonar equation can be 
written in terms of a signal-to-noise ratio ( SNR ). 

SL-TL-NL = SNR (9) 

With this arrangement detections are theoretically possible as long as the signal- 
to-noise ratio is greater than zero (e.g., signal can break out of the noise). 

The source level ( SL) is simply the sound pressure level radiated at one meter 
from the acoustic source, and the transmission loss (TL) is the reduction of that sound 
level due to propagation losses. As mentioned previously, the ray tube model describes 
transmission loss according to geometric spreading of the acoustic energy, but it is not 
the only important physical process involved. 

In addition to the geometric spreading of a sound signal, transmission loss also 
includes losses due to absorption, scattering, and diffraction. Collectively, these represent 
loss due to attenuation (Urick 1983). Absorption, also known as volume attenuation, is 
the conversion of acoustic energy into heat energy that is transferred to the seawater. It 
involves a frequency-dependent response to fluid viscosity, as well as magnesium sulfate 
and boric acid molecules found in seawater (Francois and Garrison 1982). As a result, its 
influence on transmission loss cannot be mitigated. If desired, however, the Bellhop 
model can factor in the impact of seawater absorption when performing transmission loss 
calculations. 

Losses due to scattering from a rough sea surface and diffraction are minimized 
when a receiver is within range to detect signals along a direct path or reliable acoustic 
path (RAP) and will not be covered in this discussion. However, surface reflections and 
scattering are an important concern when attempting to optimize acoustic modem 
performance in shallow water environments. 

Finally, if any excess acoustic energy remains after the impact of transmission 
loss an acoustic receiver still must contend with the noise level to make detections. Noise 
is the combination of background, or ambient, environmental noise and receiver self- 
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noise. Underwater ambient noise is commonly generated by surface winds, precipitation, 
shipping vessel density, biologies, and seismic activity (Urick 1983). It can be filtered out 
at a receiver if it does not impact the portion of the frequency spectrum that contains 
the signal of interest. Self-noise, however, can be a greater challenge especially for 
underwater gliders that were not originally designed as passive acoustic platforms. 
Nevertheless, the predominantly silent operation and mobility of gliders make the 
vehicles very appealing for acoustic sensing research. 
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IV. ACOUSTIC SENSING EXPERIMENT 


An experiment was conducted in February 2015 on the north side of Monterey 
Canyon at Smooth Ridge to assess the acoustic recording performance of an Unmanned 
Underwater Vehicle (UUV) as it transitioned into a deep loiter. A medium frequency 
underwater acoustic transducer provided the source signals and a Spray glider, equipped 
with an acoustic recorder, served as the mobile sensing platform. Previous research 
(Baumgartner et al. 2014) demonstrates that autonomous underwater gliders can detect 
and classify underwater marine mammal calls when given the capability to process 
acoustic signals onboard. The goal in this experiment was to intentionally broadcast 
continuous wave (CW) signals and evaluate their reception throughout the descent, deep 
loiter, and ascent dive cycle phases of an autonomous underwater glider. Experiment 
geometry was planned to demonstrate that a mobile glider can relocate from the acoustic 
shadow zone in order to exploit reliable acoustic propagation paths. To simplify the 
design requirements for this demonstration, onboard processing capabilities were not 
implemented. All recorded data was processed following the recovery of the glider. 

A. THE BELLHOP ACOUSTIC PROPAGATION MODEL 

The Bellhop model was used during this study to produce acoustic propagation 
predictions that supported planning the experiment geometry, comparing observed 
transmission loss with calculated theoretical levels, and assessing acoustic 
communications performance. The model is designed to address the previously described 
fundamental physics problems of underwater acoustics, namely, what path will sound 
take underwater and what will happen to the original acoustic energy along the way. 

Formally, Bellhop is a Gaussian beam tracing model that was developed to solve 
acoustic wave propagation problems and predict acoustic pressure fields in the ocean 
(Porter and Bucker 1987). For a given frequency, it uses a fan of beams to approximate 
an acoustic source and traces those beams as they propagate throughout a simulated 
underwater environment. The beams themselves are the model for how the original 
energy, or acoustic pressure, decays as it propagates through the ocean. This is 
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accomplished by associating each acoustic ray path with a beam model and a 
corresponding energy profile that decays according to a Gaussian distribution (Porter and 
Bucker 1987). This approach enables the calculation of acoustic pressure values for a 
range of receiver locations by “summing the contributions of each of the individual 
beams” (Porter and Bucker 1987). These design features allow the Bellhop model to 
produce better predictions than basic ray tracing methods and still be computationally 
efficient (Porter 2011). 

1. Model Inputs 

Input data is provided to the model in an environmental file (see Appendix A) that 
contains, at minimum, a sound speed profile and descriptions of the ocean surface and 
bottom boundaries. This file also allows a user to set the acoustic source frequency, 
enable the calculation of frequency-dependent volume attenuation, and define source and 
receiver geometry. When deemed necessary the model has the additional options to 
implement range-dependent sound speed and bottom profiles, as well as define ocean 
surface roughness. 

a. Sound Speed Profile 

The Spray glider used in this experiment is equipped with a pumped conductivity, 
temperature and depth (CTD) sensor that records fine resolution data during the ascent 
portion of its dives. After the completion of a dive, these observations were used to create 
the sound speed profile (SSP) for use in the Bellhop model environmental input file. This 
was accomplished through the use of a MATLAB script (Morgan and Pender 2003) that 
converts the temperature, salinity, and pressure measurements to sound speeds using the 
UNESCO “Sound Speed in Seawater” polynomial (Fofonoff and Millard 1983). 

An initial SSP was created with data from a dive conducted the day prior to the 

experiment. It was used to generate the acoustic transmission loss predictions necessary 

for selecting an ideal source-to-receiver range to meet our objectives. Subsequent Bellhop 

model mns were conducted using the SSP created from data collected during the 

experiment dive. This was done to simulate an acoustic propagation model that could 

theoretically run aboard an autonomous vehicle and perform calculations with real-time 
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oceanographic data. All model runs were conducted with range-independent sound speed 
profiles. This approach uses a single sound speed profile to model the ocean between an 
acoustic source and receiver. While inappropriate for long range problems, this was 
considered an acceptable approximation given the small scale of the experiment ( < 6 km 
propagation paths) and the relatively level ocean bottom. 


Figure 8. Bellhop input Sound Speed Profile (SSP) 

Sound Speed (m/s) 



The solid line indicates the SSP derived from Spray glider observations the day prior to 
the experiment after the resolution was smoothed to 10 m using linear interpolation. The 
dashed line is a February SSP from the GDEM database that was used to extend the 
glider profile to the local bottom depth. 


Testing revealed that a 10 m sub-sampled version of this SSP using linear 
interpolation provided adequate model resolution and eliminated undesirable and 
improbable artifacts in the model predictions (e.g., multiple fine scale ducts and 
channels). Smoothing of a SSP in this way is recommended in the Bellhop literature as a 
means to mitigate the influence of small localized SSP features on output predictions as 
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well as manage model run time (Porter 2011). Additionally, glider dives achieved 
maximum depths short of the bottom so the sound speed profiles were extended with 
climatological data from the U.S. Navy Generalized Digital Environmental Model 
(GDEM) database maintained by the Naval Oceanographic Office. 

b. Ocean Surface and Bottom Boundaries 

Acoustic energy interactions with the ocean surface and bottom induce signal 
direction changes according to the reflectivity of the boundaries. Given a particular 
source-to-receiver geometry and environmental conditions these boundary interactions 
can result in a complicated arrival structure. These multipath arrivals (see Figure 9) create 
interference effects at a receiver such as amplitude fading and signal distortion that can 
present challenges for both passive signal detection and acoustic communications. The 
Bellhop model can more accurately predict the presence and parameters of multipath 
arrivals when provided with realistic or range-dependent bathymetry data and ocean 
surface shape estimations. 

To tailor the model predictions in this study to the local environment the initial 
and terminal Spray glider dive coordinates were used to extract along-track bathymetry 
from a U.S. Geological Survey (USGS) database (Maher et. al 2001). This information 
was then used in a range-dependent Bellhop bathymetry file to define the ocean bottom. 
While surface-reflected and bottom-reflected acoustic rays were no doubt present during 
the in-field test, experiment geometry was specifically chosen to demonstrate the 
transition between the shadow zone and the lower loss downward refracting direct path 
rays. Nevertheless, estimates of ocean bottom geo-acoustic parameters were incorporated 
into the Bellhop environmental input file for completeness (see Appendix A). The option 
to alter the ocean surface shape in order to approximate surface roughness, or variability, 
was not utilized in the modeling mns. Instead, the default flat representation of the ocean 
surface was used. However, a small degree of transmission loss from surface interaction 
was expected due to the forecast light wind conditions. 
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2. Assumptions and Limitations 

While the ray tracing capability of Bellhop greatly assisted with planning the 
experiment, it is necessary to address two model limitations and assumptions that were 
made. First it is important to emphasize that the output of a Bellhop model run is a 
collection, or ensemble, of deterministic calculations. This means that given the set of 
user-specified inputs, the model will produce the same output every time. However, the 
unique ensemble of deterministic calculations does not represent a deterministic forecast. 
Instead, it is an estimation of the likely behavior of acoustic energy given the user- 
provided i ni tial conditions (Porter 2011). This is because the initial conditions are 
themselves a snapshot representation of small scale ocean features that randomly change 
over time. 

In other words, because of environmental uncertainty it is difficult to achieve 
truly accurate predictions from a single run of a model such as Bellhop. To incorporate 
random variation, such as surface roughness, the model would need to be run many times 
with variable input data. A probabilistic, or stochastic, forecast could then be generated 
from a statistical analysis of the resulting ensemble of model runs. While the 
development of efficient statistical acoustic channel simulators is common, particularly in 
acoustic modems research (Llor and Malubres 2013, Dol et al. 2013), we will assume that 
individual Bellhop predictions are sufficient for our purposes in this study. 

Secondly, a Bellhop prediction is a two-dimensional representation of the water 
column between source and receiver. It does not account for shifting bearings between a 
source and receiver over a period of time. However, given the anticipated source-to- 
receiver range, experiment duration, and average glider velocity, separate predictions for 
individual bearing lines are not deemed necessary. We therefore assume that a single 
Bellhop model run is sufficiently accurate to estimate ray paths and transmission loss 
along all of the changing bearing lines between the source and receiver. 

B. EXPERIMENT GEOMETRY 

Acoustic transmission loss predictions were examined to identify an ideal glider 
range and depth pair, relative to the sound source, that would promote exposure to the 
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downward refracting sound (see Figure 9). Based on the planned loiter depth of 1000 m a 
source-to-receiver range of 5.5 km was selected. The ship with the sound source was then 
positioned relative to the glider dive cycle to achieve our objectives: observe the increase 
in signal intensity as the glider transitioned from the shadow zone to the region 
dominated by reliable acoustic propagation paths. A broadside aspect, with respect to the 
glider dive cycle, was chosen for acoustic signal transmissions (see Figure 10). This was 
done to simplify range estimation and facilitate a comparison between observed and 
predicted transmission loss as a function of depth. The broadside aspect also ensured 
the glider would be in position to take advantage of reliable acoustic paths during its 
deep loitering. 


Figure 9. Bellhop acoustic transmission loss prediction 



0 1000 2000 3000 4000 5000 

Range (m) 


Acoustic transmission loss (TL) output prediction for a 2 kHz sound source (circle) at a 
25 m depth with an eigenray plot overlay. An eigenray plot depicts the multiple ray paths 
that will arrive at a receiver given its specified range and depth from a sound source. The 
rays intersect (arrow) at the chosen glider depth of 1000 m and 5.5 km range from the 
sound source. The dashed lines depict the reliable acoustic path rays that have not 
interacted with the bottom boundary. Note the surface duct predicted as a result of the 
small scale features that the SSP in figure 8 describes at shallow depths. 


24 












Figure 10. Acoustic sensing experiment geometry 
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An overhead view of the location and associated bathymetry selected for the passive 
acoustic sensing experiment conducted on 24 February 2015. The experiment geometry, 
relative to Monterey Bay, is represented by the surface track approximation of the glider 
dive from start (green arrow) to finish (red arrow), and the relative location of the sound 
source (yellow circle). 


C. DATA COLLECTION 

Two Acousonde 3A recorders (see Figure 11) were used in this experiment; one 
monitored the medium frequency acoustic source levels, and the other served as the 
acoustic receiver aboard the Spray glider. The Acousonde is a miniature, self-contained, 
autonomous underwater acoustic recorder (Acoustimetrics 2015). The recorders were 
programmed prior to the experiment (see Appendix B), and the unit aboard the Spray 
glider utilized an external battery pack to achieve nearly a month of continuous acoustic 
recording prior to the test day. 
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Figure 11. Acousonde 3A underwater acoustic recorder 
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From Acoustimetrics 2015, http://www.acousonde.com/downloads/Acousonde3A_Brochure.pdf. 


1. Source 

Acoustic signals were broadcast using a U.S. Navy Type G34 mid-frequency 
transducer (see Figure 12) provided by the Underwater Sound Reference Division 
(USRD) of the Naval Undersea Warfare Center (NUWC). It was deployed via winch 
cable from the stern of the R/V Fulmar and powered by a shipboard amplifier. A 
MATLAB routine was used to sequence the broadcast commands to a signal generator, 
and a handheld GPS unit was placed at the stem of the vessel to establish a record of 
signal transmission locations. Source levels were monitored by an Acousonde attached 
to the winch cable. A source depth of 25 meters was chosen to simulate a shallow 
contact. The depth also serves as a reasonable proxy for an acoustic modem housed in 
Wave Glider underwater towed body. Broadcast frequencies of 2 kHz and 4 kHz were 
selected for their stable narrowband reproduction by the G34 transducer. Unfortunately, 
difficulties with the power amplifier and G34 wet mate connectors reduced the time and 
output voltage available for signal transmissions. Despite these constraints, acoustic 
signals were successfully broadcast during crucial portions of the deep glider dive 
and loiter. 
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Figure 12. Type G34 transducer ready to deploy from R/V Fulmar 



2 . Receiver 

The second Acousonde recorder was mounted inside the Spray glider science bay 

and received the signals from the G34. This specific mounting location was chosen 

because it is an extremely simple and fast method to provide an acoustic recording 

capability to a Spray glider without vehicle modification. Admittedly, this places the 

recorder in close proximity to sources of glider self-noise such as the pumped CTD and 

buoyancy engine electric pump motor. It also degrades the omnidirectional recording 

capability of the Acousonde because the science bay components, and potentially the 

pressure hull itself, can acoustically mask the hydrophone. Therefore, this is not 
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recommended as a permanent location for a hydrophone on the Spray glider platform. 
Unfortunately, unexpected R/V Fulmar repairs delayed the execution of the experiment. 
As a result, upon recovery it was discovered that the Spray glider Acousonde had reached 
its storage limit prior to the glider’s ascent back to the surface (see Figure 13). 

D. DATA PROCESSING 

An initial review of glider data revealed that the dive cycle was completed in 
6 hours and 26 minutes. It included 57 minutes of loitering at a depth of approximately 
1000 meters, and covered a surface distance of 5.8 kilometers. 


Figure 13. Spray glider acoustic experiment dive 



Tiint (GMT) 

The dashed line and circle (green) indicate when signal transmissions began after 
troubleshooting delays. The second circle (red) marks the time when the Acousonde 
aboard the Spray glider reached its internal storage capacity. 


28 




Once the glider was recovered, the data from both Acousonde recorders was 
downloaded and portions applicable to the signal broadcast times were read into 
MATLAB. The companion MATLAB script (acousonderead.m) automatically accounts 
for receiver sensitivity and amplifier gain. The script reads the recorded data, for an 
indicated period of time, and calculates acoustic pressure values in micropascals (pPa) 
(Miller 2013). Once in this format, the acoustic signal data can be transformed from the 
time domain to the frequency domain using spectral analysis. This digital signal 
processing technique uses the Fast-Fourier Transform (FFT) to estimate the average 
acoustic power of digitally-sampled data with respect to its frequency spectrum (Johnson 
1989). The relationship between acoustic power and frequency is known as power 
spectral density (PSD) and is often represented in decibels (clB) per hertz {Hz.) (Stoica and 
Moses 2005). When acoustic pressure data for a particular time interval is transformed in 
this fashion the resulting PSD time series plot format is commonly referred to as a 
spectrogram. For large data sets, a technique known as the Long-Term Spectral Average 
(LTSA) can be used to present the output of spectrogram calculations in a compressed 
format for quick identification of signals (Scripps 2015). 

With these techniques, source level estimates for the transmission frequencies 
were easily extracted from the G34 Acousonde power spectrum time series (see Figure 
14 and Figure 15). A transmission loss correction was then applied to these values 
because the Acousonde reference distance was greater than one meter from the center of 
the G34 transducer. The result is a source level time series that can be used to calculate 
observed transmission loss when compared to data recorded on the Spray glider. 
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Figure 14. Spectral analysis of source level data 




A Long-Term Spectral Average (LTSA) time series (3 second averaging window, 1 Hz 
bins, 29 kHz sampling frequency) produced from the G34 source level reference 
Acousonde. Broadcast intensity and duration of the 4 kHz and 2 kHz signals is clearly 
visible above ambient noise spectrum level (SNR > 0). 


Figure 15. G34 4 kHz signal source level 



A sample of the 4 kHz signal source level derived from the data depicted in figure 14. A 
transmission loss correction, assuming spherical spreading loss, has been applied to 
adjust Acousonde pressure data for the 1.7 meter reference distance from the G34. Note 
the stability of the source signal amplitude with a variation of much less than 1 decibel. 
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The acoustic data set from the Spray glider proved to be useful in several ways: it 
contains over three weeks of acoustic ambient noise and marine mammal recordings, as 
well as an hour and a half of received signal data from the descent and loiter phases of the 
glider experiment dive. Compared to the G34, calculating acoustic receive levels at the 
Spray glider proved to be a more involved process due to unanticipated sources of noise. 

1. Ambient Noise Filtering 

In order to verify that the Bellhop model provided good predictions for setting up 
the experiment geometry, it was necessary to compare the transmitted signal intensity to 
the intensity recorded at the receiver aboard the Spray glider. Under typical conditions we 
could simply show that the signal energy level at the receiver, for the corresponding 
broadcast frequencies, was consistently observed above the ambient noise level. 
Unfortunately, the unanticipated and powerful acoustic air gun blasts of a nearby 
geological survey vessel made this a more challenging task (see Figure 16). 


Figure 16. Geological survey noise spectrogram 



In this sample of data from the Spray glider the acoustic air gun blasts begin around the 
18 second mark and contaminate the receive level estimates in the acoustic source 
frequency band (2 kHz). 
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When fired, the broadband acoustic intensity of these air gun blasts often 
surpassed the received narrowband intensity of the signals sent by the G34. Therefore, it 
was necessary to filter out these unwanted energy spikes from the received data. A 
MATLAB script was devised to read-in data and identify times when acoustic pressures 
exceeded a threshold that corresponded to air gun blasts and echoes. The stored times, 
along with a window matched to the average noise duration, were then used to ignore 
pressure data that would have otherwise introduced an unwanted bias in transmission loss 
assessments (see Figure 17). This technique revealed that 9% of the 10-minute 4 kHz 
data set, and 20% of the 80-minute 2 kHz data set, was masked by the air gun noise. 


Figure 17. Air gun noise filtering example 



Time (s) 


Signal data (2 kHz band) from the same time period as the spectrum displayed in figure 
16 after data points that correspond to air gun impulse and echo noise have been 
identified for exclusion (red). Ambient noise is plotted as a visual representation of the 
signal-to-noise ratio ( SNR ) present at the receiver. 


Once air gun noise is identified and excluded from the data, the remaining 
ambient noise levels are quite comparable to the empirically derived predictions (Wenz 
1962) for the source frequencies, 4 kHz and 2 kHz, given the prevailing environmental 
conditions (e.g., Sea State 1). 
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Figure 18. Wenz curves for ocean noise levels 
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Typical background noise levels in the ocean as adapted from Wenz (1962) by the 
National Research Council (NRC) in Ocean Noise and Marine Mammals (2003) (after 
Discovery of Sound in the Sea 2015, http://www.dosits.org/science/soundsinthesea/ 
commonsounds/?CFID=3833629&CFTOKEN=32573164) 
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Self-Noise Considerations 


The Spray glider is typically a very quiet platform for acoustic sensing. However, 
several routine onboard operations generate significant self-noise on the recorder. Two 
mechanical sources are responsible for the majority of this broadband noise: battery pack 
adjustments and buoyancy engine operation. The glider autopilot controller, with input 
from flight control dead-reckoning algorithms, commands periodic battery pack 
movements throughout dives to correct pitch and heading errors (Sherman et al. 2001). 
Additionally, the coordinated operation of the hydraulic buoyancy pump and battery pack 
ballast movement is necessary to level off a glide descent. 

These sources of noise, common among glider designs, can have a serious impact 
on acoustic sensing. One solution, implemented aboard a Slocum glider by Rogers et al. 
(2004), is to have a quiet or “comatose” operating mode available for periods of acoustic 
data collection. In a similar fashion, the glider used in this study demonstrates the 
benefits of a silent operating mode during the deep loitering period of its dive cycle. This 
is because once established at the target depth no further pitch, roll, heading, or depth 
trim adjustments were necessary. For purpose-built acoustic gliders, careful coordination 
among acoustic and non-acoustic systems is the ideal solution if systems are to be fully- 
integrated. In this study, however, the routine self-noise that the Spray glider generates 
has been carefully identified and ignored using the previously described filtering process. 

3. Source-to-receiver Range Estimation 

Ambient noise and self-noise were examined to improve the assessment of 
observed transmission loss. In order to compare this experimentally observed assessment 
with the Bellhop model predictions it was necessary to estimate glider ranges, from the 
G34 source, throughout the dive cycle. The experiment was designed to keep the Spray 
glider at a fairly constant range from the sound source to capture the transition from the 
shadow zone into the region defined by the reliable acoustic propagation paths. 

In reality, environmental factors caused the ranges to vary throughout the 
experiment. This is due to the drift of the research vessel from wind and currents, as well 
as the underwater track the glider made between waypoints. While a GPS log aboard the 
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R/V Fulmar provides a good estimate for the sound source drift, the shape of the glider 
track must be estimated based on two reported surface positions. This is because position 
updates are not available while the glider is submerged. As previously introduced (see 
Chapter II), research demonstrates that this deficiency could be solved through 
coordination with a surface navigation reference such as a GPS-enabled unmanned 
surface vehicle (USV). Since that capability was not available in this study, the standard 
mode of Spray glider navigation will be explained. 

Before a dive, the Spray glider acquires an initial GPS position and uses a range 
and heading calculation to path plan to the next waypoint. The planning algorithm does 
not factor in anticipated drift associated with the depth-averaged ocean currents from 
previous dive cycles (Sherman et. al 2001). Instead, once the glider is below the surface, 
navigation accuracy is based on dead-reckoning that is supported with feedback from a 3- 
dimensional compass (e.g., pitch, roll, and heading) (Sherman et al. 2001). As a result, 
the GPS fixes that mark where a dive begins and ends are the only available references 
for the true dive track. While other research has been done to accurately estimate and 
optimize the underwater path a glider takes (Shim 2009), a simple interpolated straight- 
line track between the GPS surface fixes will be used here. This assumption is 
appropriate given the spatial scale of the experiment, as well as the glider’s estimate that 
the depth-averaged currents would help carry it toward its target waypoint. 

Additionally, horizontal velocities and their influence on along-track progress for 
the various phases of the dive are unknown. According to Sherman et al. (2001), Spray 
glider horizontal velocities are not measured or calculated because they depend on 
estimates of the angle of attack. As a result, the following assumptions about the dive 
geometry were made: the glider only made appreciable along-track progress when 
vertical velocities, and momentum, were highest and not while it loitered at depth; the 
glider had advanced to approximately the midpoint along the 5.8 km dive-cycle track 
when it entered its loitering phase. Given these assumptions, estimates were made of 
along-track glider positions in latitude and longitude that were then paired with the 
associated dive cycle depths and times. All depths and times which corresponded to the 
deep loitering period were given the coordinates for the dive cycle halfway point. 
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V. RESULTS 


Before the unanticipated background noise was filtered out, the raw acoustic 
pressure data from the Spray glider Acousonde was first examined for the presence of the 
transmitted signals (see Figure 19). This was accomplished using the previously 
described LTSA technique (See Chapter IV). 
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Figure 19. Spectral analysis of glider acoustic data 
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A Long-Term Spectral Average (LTSA) (3 second averaging window, 1 Hz bins, 29 kHz 
sampling frequency) produced from the Spray glider Acousonde recorder data. The 
transmitted signals are visible as horizontal traces at 4 kHz and 2 kHz. Geological survey 
air gun background noise begins during reception of the 4 kHz signal as indicated by the 
high energy vertical lines. The period of very intense acoustic energy after 20:15 is 
caused by glider self-noise associated with dive descent level-off procedures. 


The output plot clearly shows that the 4 kHz and 2 kHz signals arrived at the 
Spray glider with sufficient remaining intensity to break out of the background noise. The 
plot also reveals the relative intensity of the nearby geological survey air gun blasts, as 
well as glider self-noise associated with the descent level-off. This initial representation 
of the data provides confirmation that Bellhop model predictions supported the selection 
of suitable experiment geometry. However, in order to understand how well the model 
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works as a planning tool it was necessary to compare the acoustic transmission loss 
predictions to the levels observed during the experiment. 

A. TRANSMISSION LOSS COMPARISON 

Observed transmission loss was calculated using acoustic source levels and 
receive levels derived from the G34 and Spray glider Acousonde recorders. Background 
noise and vehicle self-noise were filtered out of the data set collected aboard the glider. 
Comparisons were then made with the Bellhop model transmission loss predictions (see 
Figure 20). To accomplish this, predicted transmission loss values were extracted from 
the Bellhop model output. The values were selected using glider depths and source-to- 
receiver range estimates that correspond to the data sample times in the observed data. 
The final comparison reveals a close agreement with the 4 kHz observations. However, 
the 2 kHz comparison reveals a bias between the model prediction and field experiment 
observations. This bias, an average of 5 decibels, is likely due to a slightly different 
acoustic pressure field structure than what is predicted by the Bellhop model. The range- 
independent sound speed profile used to generate the Bellhop model prediction is quite 
possibly the cause for this difference. Nevertheless, the comparison provides suitable 
confirmation that the model predictions were adequate for planning experiment 
geometry. 

It is interesting to note that the observed transmission loss plot for the 4 kHz 
signal has variability that occurs at approximately a 10-second period. The hypothesis is 
that the sea surface swell in the vicinity of the experiment is the cause for this periodic 
variance in transmission loss. During the experiment, the long-period swell was from the 
northwest and this resulted in wave crests that were essentially perpendicular to the 
propagating acoustic energy. Additionally, the Bellhop prediction (see Figure 9) suggests 
that direct path and surface-reflected path acoustic rays would have been the highest 
intensity arrivals at the glider location. Therefore, a sea surface with a stable swell period 
of 10 seconds could certainly have induced periodic acoustic variability in the form of 
constructive and destructive interference seen at the receiver. Ultimately, the variable 
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interference is the result of small changes in the acoustic phase and intensity of the 
primary eigenrays received at the Spray glider. 


Figure 20. Observed and predicted transmission loss 
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B. LOITER ACOUSTICS 


The glider dive cycle conducted during this experiment was planned to 
incorporate an hour-long loiter at a depth of 1000 meters. While loitering at depth is an 
advertised capability of current-generation Spray gliders, the process is not optimized for 
passive acoustic sensing. The acoustic data reveal that operation of the hydraulic 
buoyancy engine pump , as well as battery pack ballast movement, result in significant 
self-noise as a loiter is being established (see Figure 21). The selected mounting location 
for the Acousonde recorder no doubt magnifies this effect. However, once glider self¬ 
noise associated with this level-off process subsides the glider is a very quiet platform for 
passive acoustic sensing. During loitering the 2 kHz signal received at the glider 
exceeded background noise levels by an average of 12 decibels (see Figure 22). 


39 












Figure 21. Self-noise during glider level-off 
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This sample of acoustic data from the Spray glider Acousonde captures the relative 
intensity of mechanical self-noise produced as the glider established a deep loiter. Noise 
was identified and filtered from the acoustic pressure data using a MATLAB script. The 
spectrogram function was then used to transform the pressure data, and sound pressure 
levels from a 7 Hz wide band centered at 2 kHz were then extracted from the output. The 
fdtered out in-band self-noise and USGS air gun noise is plotted in red. Sound pressure 
levels from adjacent frequency bands is plotted in purple as a visual representation of 
background noise and the signal-to-noise ratio (SNR) at the loiter depth. 


Figure 22. 2 kHz signal receive level at loiter depth 
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This sample of Spray glider acoustic data depicts reception of the 2 kHz signal shortly 
before the termination of the deep loiter. Analysis was conducted in the same manner as 
is described in Figure 21. Acoustic air gun noise has been filtered out and is not shown. 
Note the in-band receive level of the transmitted 2 kHz signal compared to the ambient 
noise from adjacent frequency bands. At the 60-second mark the G34 transducer lost 
power from the shipboard amplifier. The in-band receive level subsequently dropped 
back down to magnitudes equivalent to the ambient noise. 
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C. LOITER STABILITY 

The Acousonde data reveals that gliders loitering at depth have great potential 
as persistent acoustic sensing platforms. Separate three-dimensional compass data 
from the glider dive log also suggests that deep loitering can be very stable. The entire 
level-off process, from buoyancy pump initiation to approximate equilibrium, required 
3.5 minutes. Mechanical processes required in order to establish neutral buoyancy and 
re-trim pitch and roll account for motion during the first 2 minutes. Vehicle pitch and roll 
then stabilized after an additional 1.5 minutes. It is useful to mention here that glider 
pitch attitude settled from 16 degrees nose down in the descent to 23 degrees nose down 
where it remained throughout the loiter. This is likely due to the loss of lift over the 
glider’s wings as vehicle momentum decreased. Additionally, once the glider had 
stabilized it slowly drifted 30 meters below the target depth over the duration of the 
57 minute loiter. 

While vehicle pitch and roll remained fairly constant throughout loitering the 
vehicle did experience some heading drift. The three-dimensional compass heading 
information indicates that ocean currents at depth were sufficient to induce a heading 
drift of approximately one degree per minute. A subsequent glider dive on a reciprocal 
heading was also examined for comparison. That dive also included a deep loitering 
period with similar vehicle stability and a heading drift of approximately 6 degrees per 
minute. These low heading drift values are good indications that a deep-loitering UUV 
would be a stable platform for both passive acoustic sensing and acoustic 
communications. 
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VI. UNDERWATER ACOUSTIC COMMUNICATIONS 


In an unclassified April 2015 presentation at Naval Postgraduate School U.S. 
Navy Captain Doulas G. Perry, Commander of Submarine Development Squadron FIVE 
(COMSUBDEVRON 5), presented a list of challenges and needs associated with the 
integration of unmanned undersea systems. Among the action items on the list he 
identified two of the fundamental technical limitations associated with UUV operations: 
high communications latency and data exfiltration alternatives to eliminate vehicle 
recovery (Perry 2015). These two factors are linked because the faster UUVs can 
communicate useful data to operational commanders (see Figure 23), the longer 
autonomous vehicles can persist in the environment where their sensors are needed. 


Figure 23. Spray glider with acoustic modem 



Installation of an acoustic modem transducer (black cylinder) on the tail section of a 
Spray glider UUV (from Send, U., L. Regier, and B. Jones, 2013: Use of Underwater 
Gliders for Acoustic Data Retrieval from Subsurface Oceanographic Instruments and 
Bidirectional Communication in the Deep Ocean. J. Atmos. Oceanic Teclinol., 30 , 984- 
998, DOI: 10.1175/JTECH-D-l 1-00169.1). 


In the theoretical network suggested in this study, communication through the 
ocean surface boundary requires a dedicated gateway USV. The undersea collaboration 
among distributed autonomous platforms, however, requires a reliable wireless method to 
connect the network. Established communication methods that utilize the radiofrequency 
(RF) portion of the electromagnetic spectrum work well in the Earth's atmosphere, 
but they do not produce adequate results underwater. This is because, compared to sound, 
RF energy attenuates rapidly in seawater (Urick 1983). As a result, underwater 

communications using acoustic modems has become the accepted method to achieve 
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connectivity at practical ranges and depths. However, it is not a perfect solution. The 
previously discussed factors that impact the propagation of acoustic energy underwater 
such as sound speed, geometric spreading, and attenuation (see Chapter III) can 
significantly limit the speed and reliability of acoustic communications. Therefore, it is 
important to understand how underwater acoustic communications works, and why it can 
sometimes fail to achieve adequate connectivity. 

A. ACOUSTIC DATA TRANSMISSION 

An underwater acoustic telemetry modem transmits digital information, or data, 
through water to a receiver modem using acoustic signals. The transmitting modem 
encodes, or modulates, a data message to allow for underwater acoustic transmission. The 
receiving modem then decodes, or demodulates, the incoming signal to reconstruct the 
original data message (Nortek International 2015). An individual data message is a 
sequence of symbols, and each symbol is constructed with multiple bits of information 
(e.g., a binary value of either zero or one). When transmitted, symbols have a particular 
duration and this determines the number of symbols that are sent per second. Modem 
output is commonly represented by the rate of data transmission. The number of bits 
transmitted per second is known as the modem bit rate, and the number of different 
symbols transmitted per second is referred to as the modem baud rate. Ultimately the 
performance of the acoustic link between two modems depends on more than just the 
output of the transmitting modem. The physical behavior of sound in the underwater 
acoustic channel can limit the rate that symbols are successfully received, known as 
network throughput. 

B. NETWORK THROUGHPUT 

The successful reception of underwater acoustic data at a receiving modem is a 
function of several factors that describe the physical layer of the network. First, 
transmitted data symbols must travel the required path and distance through the ocean to 
the receiving modem. Secondly, the symbols must arrive with sufficient signal intensity 
for recognition above the background noise level. Finally, the receiving modem must be 
able to demodulate, or decode, the symbols with no errors. 
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1 . 


Path and Distance 


When modulated digital data travels a distance, along a particular network path, to 
a receiving modem the time delay is known as latency. When data is transmitted in the 
form of underwater acoustic modem signals, both the acoustic propagation path and the 
latency are influenced by the change in the physical properties of the ocean along the way 
(e.g., variations in sound speed). Additionally, the latency is also a function of modem 
signal encoding and decoding delays. Variability of key physical parameters in the ocean 
significantly impacts acoustic propagation ranges and affects acoustic travel times. As a 
result, the latency of an underwater acoustic communications link is quite high by 
comparison to traditional RF or wired communications networks. Exploitation of 
underwater direct acoustic propagation paths can minimize travel times and attenuation 
due to the shorter path distances. However, modem geometry alone cannot completely 
compensate for the underlying physical limitations of the medium. 

2. Signal intensity 

When a modem communicates data across an underwater acoustic channel, the 
signal intensity at the receiving modem is primarily a function of geometric spreading 
and attenuation. This means that signal intensity decays as communication path lengths 
increase. As a result, modems must have an ideal separation with respect to range and 
depth to ensure signals can arrive with sufficient intensity to be decoded. Transmission 
loss plots, introduced in Chapter IV, are a useful tool to predict the proper geometry 
that will promote a good signal-to-noise ratio (SNR) at the receiving modem. These plots 
can be used to estimate the regions of minimal transmission loss for particular 
acoustic communications frequencies. This is an important consideration to optimize 
communications because the attenuation of acoustic signals is frequency-dependent: high 
frequencies can support high data rates but only at shorter ranges, and long range 
communications are only possible with lower frequencies. This relationship between data 
rate and frequency is why underwater acoustic communications is commonly referred to 
as bandwidth-limited. In other words, the bandwidth of acoustic frequencies available for 
data transmission decreases as range increases. Unfortunately, estimating acoustic 
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modem performance is more complicated than modeling transmission loss. Reliable 
underwater acoustic communications also depends on the ability to mitigate the influence 
of multipath interference. 

3. Minimizing error 

Given a particular undersea environment, there are many possible combinations 
of depths and ranges to achieve acoustic communications connectivity between two 
modems. Optimal connectivity is achieved when a modem can receive, decode, and 
reconstruct transmitted data symbols with minimal error. In the proposed network, one 
way to promote this condition is to use a submerged towed body (see Figure 24) to 
increase the depth of the modem connected to the USV. This will minimize the impact of 
wind-generated noise, and Doppler motion associated with the moving sea surface. It will 
also improve connectivity with deep-loitering UUVs that already benefit from an 
environment of low ambient noise. 


Figure 24. Wave Glider SV2 with towed body 



A Liquid Robotics Wave Glider with an underwater towed body that can be outfitted 
with an acoustic modem (from Integrated Wave Glider-Echosounder. NOAA Fisheries: 
Northwest Fisheries Science Center, http://www.nwfsc.noaa.gov/news/ 
features/wave_glider/index.cfm). 
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However, the primary source of error in underwater acoustic communications is 
multipath interference. Acoustic reflections from the sea surface and ocean bottom result 
in multiple time-delayed arrivals of the same data symbol. If a modem cannot adapt to the 
particular time spread between these arrivals then it will not be able to reconstruct an 
error-free symbol (Send et al. 2013). This condition, known as Intersymbol Interference 
(ISI), is the reason why it can be challenging to establish reliable acoustic 
communications in high multipath environments with closely-spaced acoustic arrivals 
(e.g., shallow water). As a result, the receiving modem will need to request 
retransmissions of the data. This compensation for high bit error rates is costly in terms of 
both the time and additional power needed for acoustic transmissions. The degrading 
effects of ISI are best mitigated in two ways: exploiting direct acoustic propagation paths, 
and modem signal processing designed to adapt to local conditions. 

C. CHANNEL IMPULSE RESPONSE 

A common signal processing technique used to improve modem performance 
involves measuring the impulse response of the acoustic channel between two modems. 
In signal processing the impulse response of a system is its response to a unit-impulse or 
brief signal (Johnson 1989). When a modem transmits acoustic signals underwater it 
results in a sequence of time-delayed arrivals at the receiving modem because of the 
multipath dispersion of the acoustic energy. The arrivals also vary in amplitude and phase 
due to ocean surface and bottom reflections. The combined effect at a receiver is time- 
varying constructive and destructive interference. Given these challenges, acoustic 
modem manufacturers have developed several techniques to promote successful 
connectivity. 

For example, the Woods Hole Oceanographic Institution (WHOI) Micromodem-2 
uses an adaptive decision feedback equalizer (DFE) to automatically respond to 
variations in the underwater acoustic channel (Gallimore et al. 2010). The equalizer uses 
received signals to estimate the channel impulse response, and then calculates the 
equalizer filter weights that will optimize the communication link (Preisig 2005). In 
addition to adaptive features, manufacturers commonly use forward error correction 
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(FEC) codes to improve performance (Evologics 2015, Teledyne Benthos 2015). This 
technique uses specially coded redundant data transmissions that enable a receiving 
modem to correct errors in received data. As a result, a receiving modem only needs to 
request retransmission of received data that it is unable to correct. Regardless of the 
particular technique, all active efforts to mitigate multipath interference reduce the 
effective bit rate. Therefore, increased acoustic communications reliability comes at the 
cost of decreased network throughput (Nortek International 2015). 

D. EXPERIMENTAL RESULTS 

Building a fully-functioning version of the proposed USV/UUV network was not 
feasible given the constraints of this study. However, it was still vital to demonstrate the 
previously described acoustic communication concepts in a representative oceanographic 
environment. To accomplish this, an experiment was conducted to simulate underwater 
acoustic communications among distributed autonomous platforms. Two acoustic 
modems were used in the test: one was positioned at various depths to represent a mobile 
UUV, and a second shallow modem (6 m depth) was deployed at various ranges from the 
first to simulate a USV gateway platform. A test site in Monterey Bay was chosen with 
similar bathymetry to the passive acoustic sensing experiment conducted earlier in the 
year (see Figure 25). 

Once again, the Bellhop model and a GDEM-derived sound speed profile were 
used to produce acoustic propagation predictions. However, this time the incoherent 
representation of the acoustic pressure field was used to assess the behavior of 
propagating underwater sound. The incoherent option, when selected, ignores the 
influence of acoustic phase in order to present a more averaged representation of 
transmission loss (see Figure 26). This option is appropriate because the precise nature of 
the acoustic pressure field interference pattern is not useful given modem modulation 
schemes and the multipath time spread of arrivals (Porter 2011). Instead a simple 
representation of signal intensity with respect to range and depth is sufficient. 
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Figure 25. Acoustic communications experiment geometry 
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An overhead view of the location and associated bathymetry selected for the acoustic 
communications experiment conducted on 27 August 2015. 
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Figure 26. Acoustic modem transmission loss prediction 
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An incoherent transmission loss prediction for an acoustic modem signal centered at 14 
kHz. The plot is useful for identification of regions where signal intensity should promote 
reception. However, even with sufficient signal intensity above the background noise 
modems must also be able to mitigate the influence of multipath interference. 

Observations from this experiment, at several combinations of ranges and depths, 
resulted in connectivity between the modems. When the modems were able to 
successfully communicate, their geometry was consistent with the regions of lower 
predicted transmission loss. As expected, communication was ineffective when one 
modem was located in the region of the plot that depicts the acoustic shadow zone. In 
several instances, when connectivity issues were believed to be the result of multipath 
interference, a change in depth separation between the modems was enough to improve 
communications. The estimated acoustic ray paths and channel impulse response for two 
representative scenarios from the field test are presented in the following figures (see 
Figures 27, 28, 29, and 30). The method for plotting the Bellhop acoustic ray amplitude 
delay information is derived from Torres (2007). 
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Figure 27. Multipath arrivals for geometry that supported connectivity (source at 250 

m, receiver at 6 m) 
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A Bellhop eigenray plot that depicts estimated acoustic ray paths between modems for a 
geometry that resulted in excellent connectivity during the field test. The red line is the 
direct path, the black line indicates a path with a single boundary interaction (surface 
reflection), and the blue lines represent bottom-reflected ray paths. 


Figure 28. Channel Impulse response (source at 250 m, receiver at 6 m) 
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Multipath arrival time delays (horizontal axis) and dimensionless amplitudes (vertical 
axis) from a Bellhop acoustic channel impulse response prediction. The plot is an 
estimate of the arrivals seen at the receiver for the most favorable acoustic modem 
geometry tested during the experiment (source at 250 m transmitting to a receiver at 6 m). 
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Figure 29. Multipath arrivals for geometry that did not support connectivity (source 

at 50 m, receiver at 6 m) 



A Bellhop eigenray plot that depicts estimated acoustic ray paths between modems for a 
geometry that did not result in successful modem communications during the field test. 
Note the lack of direct acoustic propagation paths in the prediction. 


Figure 30. Channel impulse response (source at 50 m, receiver at 6 m) 
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Predicted acoustic multipath arrival time delays and amplitudes for a modem separation 
that did not allow for any successful communications during the field test. Note the 
amplitudes and time spread of arrivals compared to Figure 28. 
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These observations suggest that it should be possible to predict modem geometry 
that will promote optimal acoustic communications given the unique parameters of a 
particular underwater acoustic channel. Therefore a logical addition to the proposed 
USV/UUV system is to develop a real-time network connectivity model to optimize 
acoustic communications (Homer and Xie 2012). Ideally, the model would mn aboard 
the USV and utilize in-situ oceanographic observations to better tune network 
communications to the local environment. 

E. POWER REQUIREMENTS 

While the underwater environment imposes significant limitations on acoustic 
communications, the most important limiting factor for battery-operated UUVs is power. 
Acoustic communications requires power for both receive and transmit modes. However, 
the amount of power that is required to operate a modem transducer for acoustic 
transmissions is much more than is needed for a passive receive mode (Send et al. 2013). 
For example, the LinkQuest UWM3000H has an advertised maximum operating range of 
6000 meters with receive mode and transmit mode power requirements of 0.8 watts and 
12 watts, respectively (LinkQuest Inc. 2015). While this particular modem might not be 
the best choice for a glider UUV the listed specifications reinforce the relationship 
between power consumption and working range. Additionally, a quick survey of modem 
designs reveals that increased ranges also require larger acoustic transducers and 
consequently a larger modem form factor. This means that there is no one size fits all 
solution. Regardless of the modem, acoustic transmit power requirements will require a 
significant share of the power budget. While this is not a concern for a USV that can use 
solar energy to replenish its power supply, a battery-operated glider UUV must be 
carefully designed to communicate in a power-efficient manner. This constraint makes a 
deep glider loiter, that limits ballast movement and buoyancy engine pumping, an even 
more appealing aspect of the proposed network model. 

F. IMPROVING PERFORMANCE 

Despite the factors that influence the reliability of acoustic communications there 
are still many ways to encourage consistent levels performance. Broadly speaking, 
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performance improvements can be broken down into three categories: hardware and 
design, software, and message prioritization. 

1. Hardware and Design Features 

Two common examples of how to improve underwater acoustic modem 
performance are to utilize directional transducers and acoustic baffles (LinkQuest Inc. 
2015, Nortek International 2015). A directional transducer sends and receives acoustic 
signals over a limited range of bearings in order to focus signal intensity in the direction 
of a particular receiver. Directional transmissions can be especially useful for limiting 
multipath interference when the relative bearing to an intended receiving modem is 
known. Acoustic baffles are design features of modem installation that are used to 
physically block reception of acoustic arrivals from undesired directions. Additionally, a 
network design feature that has potential to improve connectivity is the use of purpose- 
built relay node UUVs. A vehicle specifically designed to loiter and relay data would 
extend communications ranges and allow distributed assets to connect with the surface 
gateway USV at greater ranges. 

2. Software 

The previously described acoustic communications experiment highlights 
the value of implementing onboard network connectivity models to optimize 
communications. The probability that the proposed network would remain connected 
ultimately depends on the optimal use of acoustic communications for given 
environmental conditions. Acoustic channel simulators that could predict both individual 
modem performance, as well as the performance of the greater network, would be a very 
useful tool to support network autonomy. 

Performance gains can also be realized through the use of separate models to plan 
sensor employment. This need for multiple models is necessary because optimal 
communications ranges and sensor coverage ranges will not always be similar. As a 
result, the network will need a computational method that can strike a balance between 
connectivity and sensor coverage. 
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3. Message Prioritization 

Acoustic modem transmissions require a large share of a battery-operated UUVs 
power budget. It is also not feasible to transmit all sensor data due to the bandwidth and 
data rate limitations imposed by the underwater acoustic channel. Therefore, acoustic 
communications should be prioritized in order to maximize vehicle endurance. Several 
message types will be required on a routine basis, but they will each have their own 
unique intervals. Some examples of necessary messages include contact reports, 
navigation updates, raw sensor data samples, C2 instructions, and vehicle health 
monitoring status updates. Given the variety of desired information, it might be very 
useful to explore ways to simplify acoustic messages and establish an optimal hierarchy 
for transmissions. 
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VII. CONCLUSION 


Experimental results from this study demonstrate that ocean glider Unmanned 
Undersea Vehicles (UUVs) are effective mobile platforms to support persistent acoustic 
sensing. Additionally, observations and background research suggest that underwater 
acoustic communications is a feasible method to connect an unmanned surface vehicle 
(USV) and a submerged glider UUV. 

During a passive acoustic sensing test a medium frequency acoustic source was 
deployed at a range of 5.5 km from a receiver on board a Spray glider UUV conducting a 
1000 m dive. Transmitted acoustic signals were successfully received while the glider 
was in motion, and reception improved during a quiet 1000 m loitering period. Acoustic 
propagation modeling supported the selection of experiment geometry that successfully 
demonstrated the value of a mobile sensing platform: gliders can adapt to local 
environmental conditions in order to exploit acoustic energy that propagates along the 
reliable acoustic path (RAP). It was also noted that deep UUV loiters are operating 
regimes that exhibit low platform self-noise and have the potential to be very stable. It is 
easy to see how this stability can be leveraged to maximize time on station and sensor 
coverage area. 

A separate acoustic communications experiment was conducted in similar 
environmental conditions to simulate communications among distributed autonomous 
platforms. Two acoustic modems were utilized in order to represent a UUV at various 
loitering depths, as well as a surface gateway USV. Results confirm that modem ranges 
and depths that promote ideal communications can be effectively predicted with in-situ 
observations of oceanographic parameters. However, in future applications the ability to 
maintain optimal modem geometry, as well as manage power requirements will be 
primary performance constraints. 
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A. RECOMMENDATIONS 

If networks of autonomous gliders are to function effectively in remote operating 
areas, with limited human oversight, they will need to have the ability to conduct critical 
signal processing tasks onboard. While post-processing of sensor data was acceptable for 
the demonstration in this study, it is not the ideal approach for the proposed network. A 
submerged vehicle will not be able to communicate a large quantity of raw sensor data 
over appreciable ranges in a timely fashion. Providing vehicles with organic signal 
detection, classification, and localization capabilities should be the top priority for future 
design efforts. Additionally, onboard signal processing routines need to utilize available 
power and internal data storage capacity efficiently in order to take advantage of the 
long-endurance of glider USV/UUV platforms. 

Research should also continue to consider ocean glider platforms for deep 
loitering. The more physically benign environment of the deep ocean is a better 
environment for passive acoustic sensing and acoustic communication to surface 
platforms. Increased depths also inhibit the bio fouling that leads to undesirable 
hydrodynamic drag and degraded performance. In order to enable this manner of 
operational use, it will be necessary to provide updated navigational information to 
submerged vehicles. We feel this is one of the more critical areas of focus for future 
research efforts. 

B. FUTURE PLATFORM CHARACTERISTICS 

In addition to onboard low-power signal processing future platforms designed to 
support the proposed network should have the following characteristics. First, vehicles 
should have onboard numerical models that use in-situ observations to improve 
prediction efforts. These models are necessary in order to optimize sensor placement and 
promote acoustic communications connectivity. Consideration should be given to running 
such models aboard the USV to take advantage of the renewable power supply. 
Naturally, given the battery-powered nature of current glider UUVs, onboard processing 
and modelling will greatly benefit from future advancements in power cell design and in¬ 
field recharging capabilities. 
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Secondly, the USV gateway platform should be able to provide navigation 
updates to each of the distributed UUVs under its control using underwater acoustic 
communications. This will allow glider UUVs to minimize trips to the surface and 
allocate more of their power budget to onboard sensors, signal processing, and data off- 
boarding via acoustic communications. However, as a backup capability, glider UUVs 
should still retain the option to conduct traditional surfacing for communication and 
navigation updates. This is necessary in order to mitigate the risks associated with a 
catastrophic loss of the network gateway USV due to a ship strike, storm, or interception. 

Finally, efforts to improve acoustic sensors should focus on reducing power, size, 
and weight requirements. Vector sensors and small volumetric arrays should be 
investigated for use aboard platforms tasked with passive acoustic sensing. Given their 
compact size, as well as the ability to resolve bearings to an acoustic source, these 
sensors would add great flexibility to a sensing network. If individual UUV platforms 
could resolve the location of an acoustic signal then it might be possible to reduce the 
required number of sensing nodes in a distributed network. 
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APPENDIX A. BELLHOP INPUT ENVIRONMENTAL FILE 


'Spray Prediction SSP ext. w/GDEM, 10 m interpolation)' ! TITLE 

2000.0 ! FREQ (Hz) 

1 ! NMEDIA 


'SVWT' 
51 0.0 

0.0 
0.8 
10.8 
20.8 

30.8 

40.8 

50.8 

60.8 

70.8 

80.8 
90.8 

100.8 
110.8 
120.8 

130.8 

140.8 

150.8 

160.8 

170.8 

180.8 

190.8 

200.8 
210.8 
220.8 

230.8 

240.8 

250.8 

260.8 

270.8 

280.8 

290.8 

300.8 

310.8 

320.8 

330.8 

340.8 

350.8 

360.8 

370.8 

380.8 

390.8 

400.8 

410.8 

420.8 

430.8 

440.8 


! (SSP spline, surf, vacuum,dB/wavelength, 
1261 ! DEPTH of bottom (m) 

1504.62 / 

1504.62 / 

1504.00 / 

1504.06 / 

1504.17 / 

1497.12 / 

1493.58 / 

1490.78 / 

1489.76 / 

1488.80 / 

1488.18 / 

1488.14 / 

1488.04 / 

1487.95 / 

1487.57 / 

1487.16 / 

1486.99 / 

1487.28 / 

1487.28 / 

1487.23 / 

1487.01 / 

1486.98 / 

1486.91 / 

1486.76 / 

1486.70 / 

1486.25 / 

1486.01 / 

1486.67 / 

1486.10 / 

1485.94 / 

1485.43 / 

1485.42 / 

1485.56 / 

1485.28 / 

1484.76 / 

1483.75 / 

1483.77 / 

1483.13 / 

1483.03 / 

1482.77 / 

1482.45 / 

1482.45 / 

1482.54 / 

1482.81 / 

1482.77 / 

1482.44 / 


Thorpe) 
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450.8 

460.8 

470.8 

480.8 

490.8 

500.8 

510.8 

520.8 

530.8 

540.8 

550.8 

560.8 

570.8 

580.8 

590.8 

600.8 
610.8 
620.8 

630.8 

640.8 

650.8 

660.8 

670.8 

680.8 

690.8 

700.8 

710.8 

720.8 

730.8 

740.8 

750.8 

760.8 

770.8 

780.8 

790.8 

800.8 
810.8 
820.8 

830.8 

840.8 

850.8 

860.8 

870.8 

880.8 

890.8 

900.8 

910.8 

920.8 

930.8 

940.8 

950.8 

960.8 

970.8 

980.8 

990.8 


1482.28 / 
1482.00 / 

1481.88 / 
1481.86 / 
1481.95 / 

1481.89 / 
1481.57 / 
1481.64 / 
1481.25 / 

1481.36 / 
1480.97 / 

1480.92 / 
1481.08 / 

1481.37 / 
1481.53 / 

1481.63 / 
1481.67 / 
1481.61 / 
1481.50 / 

1481.64 / 
1481.80 / 

1481.92 / 
1481.73 / 

1481.41 / 

1481.35 / 

1481.37 / 
1481.40 / 

1481.35 / 
1481.40 / 
1481.33 / 

1481.28 / 

1481.19 / 
1481.12 / 

1481.20 / 

1481.28 / 
1481.27 / 

1481.29 / 
1481.24 / 
1481.18 / 

1481.29 / 

1481.30 / 

1481.30 / 

1481.30 / 

1481.24 / 
1481.22 / 

1481.25 / 

1481.36 / 
1481.29 / 

1481.42 / 
1481.56 / 
1481.72 / 
1481.86 / 
1481.83 / 

1481.92 / 
1481.95 / 
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1000 

. 8 

1481. 

94 

/ 

1010 

. 8 

1482 . 

05 

/ 

1020 

. 8 

1482 . 

11 

/ 

1030 

. 8 

1482 . 

13 

/ 

1040 

. 8 

1482 . 

16 

/ 

1050 

. 8 

1482 . 

19 

/ 

1060 

. 8 

1482 . 

21 

/ 

1070 

. 8 

1482 . 

24 

/ 

1080 

. 8 

1482 . 

27 

/ 

1090 

. 8 

1482 . 

29 

/ 

1100 

. 8 

1482 . 

32 

/ 

1110 

. 8 

1482 . 

39 

/ 

1120 

. 8 

1482 . 

4 6 

/ 

1130 

. 8 

1482 . 

53 

/ 

1140 

. 8 

1482 . 

60 

/ 

1150 

. 8 

1482 . 

66 

/ 

1160 

. 8 

1482 . 

73 

/ 

1170 

. 8 

1482 . 

80 

/ 

1180 

. 8 

1482 . 

87 

/ 

1190 

. 8 

1482 . 

94 

/ 

1200 

. 8 

1483 . 

01 

/ 

1210 

. 8 

1483 . 

09 

/ 

1220 

. 8 

1483 . 

17 

/ 

1230 

. 8 

1483 . 

25 

/ 

1240 

. 8 

1483 . 

33 

/ 

1250 

. 8 

1483 . 

41 

/ 

1260 

. 8 

1483 . 

48 

/ 

1261 

. 0 

1483 . 

80 

/ 

■ A* ’ 

O 

O 




1261 

1503. 

0 0.0 

l.: 

352 


1 

25.0 / 

1600 

0.0 1500 / 

1600 
0.0 8 / 

'CB' 

0 

-90.0 90.0 / 
0.0 1500.0 8.0 


0.1458 / 

! NSD 

SD(1:NSD) (m) 

! NRD 

! RD(1:NRD) (m) 

! NR 

! R(1:NR) (km) 

! 'R/C/I/S' 

NBeams 

! ALPHA1,2 (degrees) 
! STEP (m) , ZBOX (m) , 


RBOX (km) 


Notes on Bottom parameters: the ocean bottom is modeled as an acoustic half¬ 
space and the range-dependent bathymetry file described in Chapter IV is called on the 
same line ('A*'). The subsequent parameters indicate depth of the bottom (1261 m), 
compressional sound speed of the bottom sediment (1503 m/s), bottom type density 
(1.352 gm/cm 3 ), and bottom type attenuation coefficient (0.1458 dB / A ). 
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The bottom type (Hemipelagic Terrigenous Clay) is derived from the Bottom 
Sediment Type (BST) database maintained by the U.S. Naval Oceanographic Office 
Acoustics Division as pail of the U.S. Navy Oceanographic and Atmospheric Master 
Library (OAML). The additional bottom geo-acoustic parameters are estimates based on 
empirically derived dependencies (Hamilton 1980, Hamilton and Bachman 1982). 


64 



APPENDIX B. ACOUSONDE SETTINGS 


Spray glider Acousonde 


Serial Number 

B003A044 (Spray glider) 

Channel 

1 (Low Power) 

Sampling Frequency 

29,038 kHz 

Cutoff Frequency 

9,292 kHz 

Total path gain (antialiasing filter engaged) 

+22.5 dB 

High-pass filter cutoff frequency 

22 Hz 

Sensitivity after preamplification 

-187.3 dB re 1 V/micro Pa 


G34 source monitor Acousonde 


Serial Number 

B003A011 (G34 source monitor) 

Channel 

1 (Low Power) 

Sampling Frequency 

29,038 kHz 

Cutoff Frequency 

9,292 kHz 

Total path gain (antialiasing filter engaged) 

+2.4 dB 

High-pass filter cutoff frequency 

38 Hz 

Sensitivity after preamplification 

-188.1 dB re 1 V/micro Pa 
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